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Abstract 



We develop the derivation we proposed in hep-th/0703177 of the dressing phase of 
the S-matrix in the AdS / CFT correspondence in the framework of the underlying bare 
integrable model. We elaborate the configuration of the Bethe roots describing the 
physical vacuum, which consists of a long Bethe string stretched along the imaginary 
axis and stacks distributed along the real axis. We determine the distribution of all 
Bethe roots in the thermodynamic limit. We then directly compute the scattering phase 
of the fundamental excitations over the physical vacuum and reproduce the BHL/BES 
dressing phase. 
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1. Introduction 

Integrability has been providing us with new insights into the duality between the 
Af = 4 super Yang-Mills theory and the superstring theory in AdS 5 x S* 5 . After 
the discovery of integrability in the one-loop super Yang-Mills theory [1, 2] and in 
the classical superstring theory [3], a lot of progress has been made toward the all- 
order /quantum integrability in the full theory of the planar AdS/CFT correspondence. 
As a monumental result, there has emerged a novel integrable model [4, 5] which is 
expected to describe the spectrum of the infinitely long Yang-Mills operators as well 
as that of the infinitely long quantum strings [6-8] , at arbitrary values of the 't Hooft 
coupling constant A. 

The integrable model is characterized by the dispersion relation [9] and the S-matrix 
[5, 10, 11] of the fundamental particles. The system exhibits the centrally extended 
psu(2|2) © psu(2|2) symmetry. Remarkably, the symmetry completely determines the 
dispersion relation and also the S-matrix up to an overall scalar factor [5]. Given the 
scalar factor as a function of the momenta and the coupling, one can systematically 
study the spectrum of the system by making use of powerful techniques developed for 
conventional integrable models, such as the Bethe ansatz. 

The determination of the scalar factor, or equivalently its principal part called the 
dressing phase, was one of the main outstanding problems in this field. The entire ex- 
pression of the dressing phase as the strong coupling expansion was first constructed [12] 
so that it satisfies the crossing symmetry [13] and includes the previously known first 
two terms [14-17] which reproduce the semi-classical string spectrum. Subsequently, a 
systematic way of its determination as the weak coupling expansion was presented [18]: 
The problem can be rephrased in terms of the cusp anomalous dimension [19] where 
the transcendentality principle [20], with some empirical rules, fully determines the 
dressing phase up to an overall multiplicative constant. The constant is readily sin- 
gled out by comparison with either the perturbative computation [21] or the strong 
coupling result [12]. Ultimately the weak coupling result is identified with the strong 
coupling one by a sort of analytic continuation [18, 22] and is nicely expressed in a 
closed integral formula. We call it the BHL/BES dressing phase after the authors of 
the articles [12,18]. Its properties, such as the pole structure [23] as well as the strong 
coupling limit [24-27], have been further studied. 

Despite the success in the determination, the clear understanding of the scalar 
factor was still lacking. The above procedures do not explain why the scalar factor 
should exhibit its particular structure. It is also unsatisfactory that these procedures 
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require some model-specific computation of the string/gauge theory. Although there 
are some interesting results explaining part of its structure [23,28-30], one would desire 
a comprehensive explanation. 

Let us recall here that in the field of integrable models, there are two well-known 
approaches for the computation of the S-matrices: One is called the factorized boot- 
strap program or the phenomenological computation [31], the other is called the direct 
calculation, the microscopic derivation or the Bethe ansatz technique [32-34]. 

The former approach is to compute the S-matrices as an inverse problem. In two- 
dimensional massive relativistic integrable models, two-body S-matrices of the funda- 
mental particles satisfy the unitarity, the factorizability, and the crossing symmetry. 
These conditions constrain the form of the S-matrices up to the CDD ambiguity. The 
ambiguity can be removed by some additional requirements, such as the absence of the 
poles corresponding to unphysical particles. 

The latter approach is to compute the S-matrices as a direct problem. For exam- 
ple, in the anti-ferromagnetic Heisenberg spin-chain the physical vacuum is the anti- 
ferromagnetic state rather than the ferromagnetic state. The anti-ferromagnetic state 
is realized as a nontrivial solution of the bare Bethe ansatz equations built over the 
ferromagnetic reference state. In other words, the physical vacuum is constructed by 
filling up the Dirac sea over the bare vacuum. The R-matrix describes the scattering 
of the magnons, which are the fundamental excitations over the bare vacuum. On the 
other hand, the S-matrix appears as the scattering matrix of the spinons, which are 
the fundamental excitations over the physical vacuum^] 

The above mentioned determination [12] of the scalar factor of the AdS/CFT S- 
matrix basically followed the former bootstrap program. It is natural to expect that 
one could determine the scalar factor alternatively by the latter direct computation. 
The idea of such nontrivial structure of the physical vacuum in the context of the 
AdS/CFT correspondence has been sometimes considered [35]. A concrete hint was 
observed in a computation of all-order anomalous dimensions [36]. The authors of [36] 
derived the integral equation describing the all-order anomalous dimensions of field 
strength operators and found that there appear integral kernels very similar to those 
describing the scalar factor. Such kernels are generated by the elimination of density 

1 The R-matrix (the scattering matrix of magnons) is proportional to the S-matrix and has a 
trivial scalar factor. It also satisfies the unitarity, the factorizability, but does not satisfy the crossing 
symmetry. In the latter approach one does not assume the crossing symmetry anywhere. Instead, 
the S-matrix, which describes the physical scattering, becomes crossing-symmetric automatically, even 
though one starts from the crossing-non-invariant R-matrix. 
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functions of Bethe roots at nested levels. In our previous article [37], we demonstrated 
that a certain configuration of the Bethe roots at nested levels indeed generates the 
dressing phase in the all-order Bethe equations^ 

In this article, we investigate in detail this microscopic formulation of the AdS / CFT 
S-matrix. After a brief review of the S-matrix and the all-order Bethe equations in Sec- 
tion 2, we present in Section 3 the whole configuration of the Bethe roots describing 
the physical vacuum. The configuration consists of a long Bethe string stretched along 
the imaginary axis and stacks distributed along the real axis. The former part corre- 
sponds to the configuration of a pulsating string, while the latter is analogous to the 
vacuum configuration of the Hubbard model in the attractive case. We determine the 
density distribution of the stacks. In Section 4, we subsequently compute the density 
of stacks in the presence of fundamental excitations. Using this, we directly compute 
the S-matrix as the two-body scattering matrix of the fundamental excitations over the 
physical vacuum, to find precisely the BHL/BES dressing phase [12,18]. Section 5 is 
devoted to a discussion. The derivation of the effective momentum phase of stacks is 
presented in Appendix A. 

2. S-matrix and nested Bethe ansatz equations 

Let us start our discussion with an introduction of some notations. The S-matrix is 
most concisely expressed with the help of the following parametrizations 

x ± (u) = x(u±i), x{u) = ! (l + V 1 - 4# 2 /<u 2 ) . (2.1) 

Here u is an analogue of the rapidity parameter and 

is the normalized coupling constant. In terms of these parameters, the momentum p of 
a fundamental particle is expressed as 

e ip = — . (2.3) 

x~ 

The scattering matrix appearing in the context of the AdS/CFT correspondence 
exhibits the following tensor product structure 

S(Pk,Pj) = S (pk,Pj) 2 [R{p k ,Pj) ® R(Pk,Pj)\, (2.4) 



2 After the submission of our previous article [37], there appeared a similar computation in the 
revised version of [36] . While the mechanism of the generation of the dressing phase is essentially the 
same, their formulation looks conceptually different from ours. For instance, there appear Bethe roots 
at the nested levels twice as many kinds as ours. 
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where R(pk,Pj) is the su(2|2) invariant R-matrix of size 16 x 16 and So(pk,Pj) 2 is 
the overall scalar factor. The form of the R-matrix is completely determined by the 
symmetry [5, 10]. There are some variations of the canonical form of the R-matrix, 
depending on the choice of the basis. Here we adopt the string theory basis [11] so that 
the R-matrix satisfies the ordinary Yang-Baxter algebra. 
The overall scalar factor is conventionally expressed as 

So( Pk , Pj ) 2 = k + i V.* i ^^, (2.5) 
x k ~ x j 1 - 9 l x k x j 

where 29(uk, Uj) is called the dressing phase [14]. If we regard the S-matrix ( 12. 4p as the 
scattering matrix of physical particles, the dressing phase turns out to be a nontrivial 
function. Its form was recently determined [12,18]. Let us call it BHL/BES dressing 
phase. It is expressed as 

/OO /'OO / 

dte^e-r / dfcV^e-TT (k d (2gt, 2gt') - K d (2gt', 2gt)) , 
■OO J —OO 

(2.6) 

where the Fourier transform is a skew combination of the dressing kernel 

/•OO ±11 

K d (t, = 8g 2 / dfK.it, 2gt")^—K (2gt", t') . (2.7) 
Jo e* - 1 

The constituent kernels are given by 



k (f n = tJ l {t)J Q {t')-t'.h[t)J 1 {t') k = t , J 1 (t)J (t)-tJ (t)J 1 (t') 



where J n (t) are Bessel functions of the first kind. 

The goal of the present article is to derive this BHL/BES dressing phase in the 
context of the underlying bare integrable model. In other words, we describe the 
system starting from a bare vacuum where the scattering matrix of the fundamental 
excitations has the same structure as (I2.4p - (I2.5I) but with the trivial dressing phase 

26 h ^ c (u k , Uj ) = 0. (2.9) 

In the bare description, the physical S-matrix can be computed as the scattering matrix 
of the fundamental excitations over the Fermi surface. 

Given the form of the S-matrix, one can derive a set of Bethe ansatz equations. 
Let us consider the system of iV particles in a periodic one-dimensional box of length 
L. We impose integrability of the system, namely the condition that any multi-body 
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scattering is factorized into a product of two-body scatterings described by the above 
S-matrix. For simplicity we consider the case of zero total momentum 

JV 

i=i 

The consistency conditions for the periodicity give rise to the Yang equations 

N 

e ipkL = Y[S(p k , Pj ). (2.11) 

These matrix equations are diagonalized with the help of the nested Bethe ansatz. 
Combined with the momentum condition (12.101) . the nested Bethe ansatz equations 
can be expressed as the asymptotic all-order Bethe ansatz equations [4] 

i=n— ( 2 - 12 ) 



=1 X 4J 



l = TT 1|fe ~ 2 ' j 1 TT - 1 ' 4J (213) 

J=I U ^ k ~ U ^3 ~ l l 2 7=1 1 - f/xi,k X^j ' 

1 — TT U2 > k ~ U2 ' j ~ 1 TT U2 ' k ~ M3j ' + ^ 2 TT M 2,fc - Mi,j + i/2 ^ 

1 = || ^-«2 J -+»/2 ^ X3, fc -X+ 

y ^3,fc - w 2 j - V 2 y ^s,fc - x AJ 

x tk \ TT ^ 4,fc - Ui,j + - e 2i0( - U4 ' k ' U4 >^ TT 1 g2 ^ x ^ k XlJ TT X/t ' fc ~ X3 ' j 

x 4,* / ^4,fc - «4j - i 1 - #7 < fc fj[ x+ k - ^3,i 

n 1 g I X 4,k X 7,j TT X 4,fc X 5J ^2 16) 
,- =1 1 y I x 4,fc x 7,j . =1 J^fc x 5j 

1 = TT ^5,fc ~ «6j + V 2 TT X 5,fc ~ X 4 J ^ 1? x 

1 _ FT M6,fc ~ Uqj - i -py n 6|A; ~ U 5J + i/2 -pr n 6 , fc - M 7 ,j + i/2 ^ 

jA u 6,fc - «6 j + i jJi ue,k - Us,j - i/2 jj^ u 6:k - u 7J -i/2' 
1 = ft M7 ' fc ~ U6 ' j TT" — (2 19) 
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The length L and the number of particle N are interpreted as 



L = J-K A + -{-K x + K 3 + K 5 - K 7 ), N = K 4 . 



(2.20) 



We refer to [5, 10,38,39] for the details of the derivation. 

3. Bethe root configuration of the physical vacuum 

In the bare description, physical states are characterized by solutions of the bare Bethe 
ansatz equations, that is, the simultaneous equation fl2.12p - fl2.19p with the trivial dress- 
ing phase f)2.9p . In this section we present a particular solution that should express the 
nontrivial physical vacuum state. 

3.1. General structure 

The configuration consists of the following occupation numbers of bare Bethe rootg§ 



For the vacuum state the configuration of Bethe roots must be symmetric with respect 
to the interchange of the two su(2|2) sectors: distribution of roots U\^i u 2,k, u^k is just 
the same as that of u 7) k,UQ^,u 5i k, respectively. Regarding this symmetry, we mostly 
omit to mention the former copy of roots hereafter. 

The vacuum has to be neutral with respect to the pair of su(2|2) symmetries. This 
restricts the relative numbers of the Bethe roots to be K± = K5 + K 7 = 2K%. It can 
be understood as follows: We restrict ourselves on one of the su(2|2)'s. Let [ni;n 2 ] 
denote the two su(2) charges, by which we mean the Dynkin indices with respect to 
the bosonic subalgebra su(2) © su(2) C su(2|2). A bosonic root u± creates a magnon 
with charges [1; 0]. Either of fermionic roots 115 or u 7 converts the magnon charges [1; 0] 
to [0; 1]. A bosonic root Uq flips the latter su(2) spin down, namely it converts [0; 1] 
to [0; —1]. Either w 5 or u 7 converts [0; —1] to [—1; 0]. In other words, M5 and u 7 have 
charges [— 1; 1] while u$ has charges [0; —2]. It then follows that a state with general 
excitations has charges [K4 — K 5 — K 7 ; K 5 + K 7 — 2Kq\. 

3 The occupation numbers (|3.1[) satisfy the condition K2 < K\ + K3 < K4 > K$ + K 7 > K§ 
required for the all-order Bethe ansatz equations. (This condition follows from the consistency of 
nested Bethe ansatz. See, e.g. [38].) On the other hand, they are outside the bound K\ < K 2 < 
K3 < K4 > K$ > Kq > Kj required for the one-loop Bethe ansatz equations. This means that the 
vacuum configuration is characteristic of all-order Bethe ansatz equations and becomes singular in the 
one-loop limit. 



{K x , ...,K 7 ) = (2M, M, 0, 2M, 0, M, 2M). 



(3.1) 
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The distribution among K 5 and K 7 is not determined by the neutralness of the 
vacuum, since the Bethe roots u$ and u-j originate in the same nested level of diago- 
nalization [5] and thus carry the same su(2|2) charges. In fact, the Bethe roots x 7 ^ 
are introduced by the relabeling x 7 = g 2 /x§ [4] in order to recover the seven sets of 
equations Gnty - GHty out of five [5,38]. 

Distinction of the roots u 5 and m 7 arises in connection with the su(2, 2|4) one-loop 
Bethe equations. For a general value of u, the value of x(u) has an ambiguity of the 
square root branch (12.11) . The branch of £5 and Xj are chosen so that x^,Xj approach 
M5,«7, respectively, in the one-loop limit g — > 0. In other words, we relabel £5 and x-j 
through the relation x 7 = g 2 /x 5 so that all x 5 's and x 7 's satisfy |x 5 | > g, \x 7 \ > g. The 
vacuum configuration has no £5 root, hence all the roots at the nested levels decouple 
from 114 s in the one-loop limit. The set of occupation numbers (13.11) are not allowed 
for the one-loop Bethe equations. The vacuum configuration is characteristic of the 
all-order Bethe equations. 

When we derive the dressing phase, we send both J and M to infinity. However, 
for the purpose of studying the vacuum configuration, it is convenient to take the limit 
J — > 00 first while keeping M sufficiently large but finite. We postpone taking the limit 
M — > 00 until we discuss excited states in the next section. 

In what follows we will specify the whole configuration of the bare Bethe roots. 

3.2. Configuration of the central Bethe roots 

In our previous article [37], we discussed mostly the configuration of Bethe roots other 
than the U4 roots, which is the most essential part of the derivation of the dressing 
phase. Here, we further specify the precise configuration of the roots of the vacuum 
state. 

The configuration of u^s is extremely simple when we see it on the w-plane. For 
J — > 00, it is given by 

U4,k = ki (3.2) 
in terms of a shifted index k = k — M — |, which runs over 

k = -M + f, -M +§,..., M-|,M-|. (3.3) 

This configuration looks like nothing but a conventional Bethe string of length 2M. In 
the present case, however, it is not enough to specify only the values of w^'s because 
for each x^%. there is a choice of two branches of the square root (12. ip . We choose them 
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in such a way that Imi^ > 0, Imx 4fc < for all roots. In other words, the vacuum 
configuration is completely specified on the x-plane. Explicitly, it is given by 



X 4,k 



\{k±\±\J{k±\Y + ^. (3.4) 



Note that the distribution ( 13. 21) on the w-plane is common to the magnon bound state 
[40]. For that state, however, the choice of branches is x\ k = x^ k+1 for k = 1, ... , 2M— 1 
and Imx^ 1 < 0, lmx^ 2M > 0, which is different from (13.41) . 

Several comments are in order. First, it is very natural that the configuration is 
simple and, in particular, does not have any continuous modulus parameter. There is 
only one discrete parameter, the number of roots 2M, which will be eventually sent 
to infinity. 

Second, the vacuum configuration (13.21) for large M is transparent when scattered 
with extra roots. More precisely, when one scatters an extra root with the vacuum 
configuration (13.21) . it gains a scattering phase against each constituent u^. However, 
there occurs cancellation and thus the total scattering phase is 

n~ !/ u 4 - w 4 j + i _ u 4 + (M + \)i Ui + (M - \)i 
x u A - u 4J - i ~ Ui - (M+\)i u 4 -(M- \)i ' 

which becomes trivial in the large M limit. This property is common to the magnon 
bound states and is crucial later in the computation of the dressing phase where we in 
fact add extra 114 roots to the vacuum. On the other hand, in contrast to the case of 
the magnon bound state, there occurs no cancellation in the parts where x^s appear 
explicitly in the Bethe equations. This is necessary for having a sufficient number of 
stack solutions; otherwise such a cancellation in (I2.19p decreases the number of solutions 
of Uj t k satisfying \xt^\ > 9 less than 2M. 

Third, the vacuum solution is characteristic of the all-order Bethe ansatz equations: 
If we take the one-loop limit g — > 0, the configuration H3.4[) becomes singular, which is 
in accord with the fact that the dressing phase vanishes in this limit. This is again in 
contrast to the magnon bound state, which survives for g — >• with each pair of x\ k 
approaching M 4 fc ±i/2. 

Some readers might wonder whether the above Bethe string with the present branch 
choice really exists, though it solves the Bethe equations in the limit J —>■ 00, and how 
to understand such a singular behavior in the one-loop limit. To answer this question, 
it would be instructive to consider the configuration temporarily in the physical Bethe 
equations, where we can make use of the correspondence with classical strings. The 
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deviation due to the presence of the dressing phase is within the error of the string hy- 
pothesis and is negligible for J — > oo. Let us consider the thermodynamic limit J — > oo, 
keeping g,M proportional to J, and introduce a rescaled spectral parameter x — x/g. 
On the x-plane the imaginary roots (13 .4p form two condensates [—ib, — 
with b = M/2g + ■>/ (M/2g) 2 + 1. Because the configuration is symmetric under the 
interchange x <-> 1/5, the corresponding classical string lives in the S 2 xf sector [41]. 
There are few candidates for the solution in the ^xR sector with only one modulus 
b. We identify it as a pulsating string [42-44] (see also [45,46]). Pulsating string is an 
elliptic solution and has a continuous elliptic modulus k and a discrete winding num- 
ber. The winding number is read from the density of the imaginary Bethe roots on the 
li-plane, which is 1 in this case. Given the winding number, the elliptic modulus k is 
determined by b. When we send M to infinity, b also goes to infinity and k approaches 
oE Thus the configuration ( 13. 4ft with large M corresponds to the rational limit of the 
pulsating string. It sweeps the S 2 at almost constant speed with high frequency. 

An unusual feature of this configuration is that the condensates run across the 
unit circle on the x-plane. Such a solution is precisely an exception to the general 
correspondence between classical strings and solutions of the one-loop Bethe equations 
[42,47,48]. 

Pulsating string solution has zero angular momentum in the S 2 . This is akin to the 
neutralness of the anti-ferromagnetic state in a spin-chain. 

3.3. Formation of stacks 

In our vacuum configuration, the Bethe roots u§ and u-j form stacks [37] 

«7,afc-i = u e,k + f , u 7j2 k = ^6,fc - | , for k = 1, . . . , M. (3.6) 

Without knowing the bare Hamiltonian, one cannot verify that this configuration really 
corresponds to the ground state. However, most likely it does, by analogy with the 
Hubbard model. The Bethe equations (I2.18p - (l2.19p resemble very much the Lieb-Wu 
equations for the one- dimensional Hubbard model. The vacuum of the Hubbard model 
was well studied [49,50]. In the attractive case, the vacuum consists of precisely this 
kind of stacks [50], namely a kind of k-A strings [51]. Note that this kind of stack also 
appears in the description of the field strength operators Tr T L [36]. 

4 Although one can take k arbitrarily small, it cannot be strictly zero as far as the winding number 
is nonzero. The strictly rational case k = corresponds to zero winding number, which looks no longer 
a pulsating string but rather a point-like string. 



9 



Multiplying Bethe equations for Uj^k-i and for u-j^k together, one obtains the fol- 
lowing set of Bethe equations 



jLi 1 - #7 x tk x lj 1 - 7 ^ »4j ^ M e,fe - «6j - i " 
They can be viewed as effective Bethe equations for denoting the centers of stacks. 

3.4. Distribution of stacks 

Given the configuration of x 4k , ( 13.71) can be viewed as Bethe equations for a single kind 
of Bethe roots with a regular form of self-interaction: 



e 



M 

i*(«6.») _ TT M e,fc - u 6 ,j + 1 



where 



*(«^) = 1 e in l "lify 1 (3.9) 

z ~t 1 - / < fc <; 1 - 9 2 / x 6 , k x tj 

is regarded as the virtual momentum phase. For sufficiently large M, one can evaluate 
this phase function by approximating sum by integral. We relegate the detail of calcu- 
lation to Appendix A. If we take M and u sufficiently large compared to the coupling 
constant g, the phase function approaches a reasonably simple form 



$(u) = 2M 



u u ( M 2 

2 arctan — + — In 1 H — 

M M V u 



(3.10) 



An important property of the function <&(«) is that it is a monotonically increasing 
function. This is clear from the form of its derivative 

&(u) = 2 In (l + ^Pj- (3-11) 

The Bethe equations ( 13.81) are thus analogous to those of the one- dimensional Bose gas 
with repulsive 5-function interaction or those of the sl(2,R) spin-chain. 
The Bethe equations (13.81) can be written in the logarithmic form 

M 

2nn k = + 2 arctan(u fe — Uj). (3-12) 
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We often abbreviate u§^ as u k hereafter. The mode number n& associated with the root 
Uk takes integer/half-integer value, depending on M is odd/even, respectively. Since 
the r.h.s. is monotonically increasing as a function of u^, it follows that > rij for 
Uk > Uj. For the vacuum configuration, we consider consecutive set of mode numbers. 
One can always relabel the Uk roots so that u k > Uj for k > j. The mode numbers for 
the vacuum configuration are then given by 

n fc = -^,-^,...,^,^ for k = l,...,M. (3.13) 

In contrast to the case of the anti-ferromagnetic vacuum of the Heisenberg chain, 
the present configuration does not correspond to the maximal filling over the real axis. 
In other words, the support of the distribution of u^s is a finite interval. One can see 
this as follows: If the real axis were occupied by u%, . . . ,um, an extra real root with 
mode number would have to sit at u = oo [33]. However, for u^ = oo the r.h.s. 
of (I3.12p would take (3M — 1)tt and thus Uj < oo for 2nrij < (3M — 1)tt, which is 
contradictory to the last argument. 

We are interested in the distribution of u^s in the large M limit. From the form 
of the potential (13.101) . we see that the characteristic length of the distribution of u^s 
is of order M. Regarding this, let us expand the summand of the interaction term in 
fl3TT2|) as 

arctan(w fc - uj) = ^sign(w fc - Uj) h O ( — rr J (3-14) 

2 u k -Uj \KUk-UjYJ 

and evaluate the summation term by term. One finds that the sum of the first term 
precisely gives rise to the mode number 

1 M 

7T sigri(^fc - U j) = n k , (3.15) 

while the sum of the lower order terms after the second one becomes negligible. The 
Bethe equations (13.121) then reduce to 

M 

$K) = 2V • (3.16) 

In the continuous limit, one can replace the sum by the principal- value integral and 
obtains 

$(.) = 2-T^, (3.17) 

J-B U ~ V 
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where we introduce the density function as 

M 



(3.18) 



3=1 



As we mentioned above the density has a finite support, which is denoted by [-B, B\. 
The integral equation (I3.17P can be solved by the inverse Hilbert transformation 

" B IB 2 - u 2 <5>{v)dv 
i_ B V B 2 — v 2 v — u 



p(u) 



2tt 



(3.19) 



The endpoints ±B are determined by the normalization condition 

-B 

p{u)du = M. 



(3.20) 



—B 



By evaluating these integral expressions, one obtains the distribution of the stacks. 

For our later purpose, let us estimate the order of B with respect to M. It is 
convenient to rewrite the above equations in terms of the rescaled variables 



u = Mu, B 
We also introduce the normalized functions 



MB. 



2 arctan u + u In ( 1 + — 



u- 



M 
3=1 



u — u 



J J ■ 



which are related to the original functions by 

= M$(fi), p(u) = Mp(u). 
In terms of these rescaled quantities, the integral equation (13. 17j) gives rise to 

rB ~f~\ 



p{v)dv 

B U-V 



(3.21) 

(3.22) 
(3.23) 

(3.24) 
(3.25) 



This equation is formally the same as (13.171) . thus the solution is given by (13.191) with 
all quantities replaced by the rescaled ones. Note that M-dependence now enters only 
through the endpoint value B, which is determined by the normalization condition 

.B ! 

(3.26) 
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It is clear that B becomes small as one sends M large. This means that for large M, 
p(u) is determined by the form of <&(u) only at small u. Except for the very vicinity of 
the origin, at small u roughly behaves as a linear function 

$(u)~4ulni, (3.27) 

where Bq ~ B is a typical scale. With this approximation one can analytically solve 
the integral equation and obtains 



* V BoJ 



p(u) ~ - In ^ V B 2 - u 2 . (3.28) 



The normalization condition (I3.26j) now reads 

S 2 lni~— . (3.29) 

Ignoring the correction coming from the logarithm, one finds that B roughly scales 
with M~ x l 2 . Thus the original B roughly scales with M 1//2 . 

4. Excited states and computation of the scattering phase 

In this section we consider excited states and compute the two-body S-matrix of fun- 
damental excitations. By making use of the underlying symmetry, the S-matrix can 
be written most generally in the form of the spectral decomposition. As the symmetry 
fixes the form of the projectors, it is enough to determine the eigenvalues in front of 
the projectors. For example, the su(2) Zamolodchikovs' S-matrix is constructed by 
computing the scattering phases with respect to the triplet and the singlet. In the 
present case, the centrally extended su(2|2) algebra possesses a peculiar feature that 
the tensor product of a pair of 4-dimensional atypical representations is irreducible [10]. 

y one scattering phase of a pair of fundamental 
This allows us to restrict our consideration of 



Therefore it is enough to compute on 



excitations in a representative stated 
excited states to those with only ii 4 roots added. 

5 The computation presumes that the physical vacuum is a singlet. In the last section we con- 
structed the vacuum as a neutral state under the pair of su(2|2) symmetries. For the centrally ex- 
tended algebra, however, the state has to be neutral with respect to the central charges as well. We 
define the action of the central charges in our bare integrable model so that the physical vacuum has 
zero central charges, by shifting one of the central charges of the reference vacuum. 
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4.1. Fundamental excitations 



Let us consider excited states by adding extra w 4 roots to the vacuum configuration. 
We let w^k denote the extra roots and A/4 be their total number. The occupation 
numbers read 

(K u ...,K 7 ) = (2M, M, 0, 2M + N 4 , 0, M, 2M). (4.1) 

We keep the structure of the other roots unchanged, namely the other 2M m 4 's con- 
stitute the Bethe string with the branch choice (13 ,4p . and the UqS and the Uj§ form 
M stacks with consecutive mode numbers. In this subsection let us determine the 
deviation of the density of stacks for fixed w 4 's. 

The effective Bethe equations for the centers of stacks read 

M 

2itn k = $ cx (wfc) + 2 arctan(w fc - Uj), (4.2) 

with mode numbers (I3.13p . The only difference from (13.121) is the momentum phase 

= $(«) + (p(u), (4.3) 
where the modification part is given by 

<^Kfc) = -> in- — ■ ; J - — ■ ; J (4.4) 

with yfj = x ± (w4j). By taking the derivative with respect to Uk, (14.21) gives rise to 

2vrp ex H = KM + 2 r 1 PcM J V ; , (4-5) 



Be: 



(u-vY + l 



where p ex (w) is the density function of the centers of stacks under the modified momen- 
tum phase. Apparently, one could derive the same integral equation for the vacuum 
density p(u) with the vacuum phase $(«). Subtracting it from (14.51) . one obtains an 
integral equation for the density deviation a(u) = p ex (u) — p(u), as follows 

f B <r{v)dv ( r B [ Bc *\ p cx (v)dv , A , 

We consider the large M limit keeping the other parameters u> 4j fc, A/4 fixed. In this 
case the second integral is negligible: The fluctuation of the overall shape of $(«) should 
be suppressed at most within the change of M to M + AM with AM = 0(1). Then 
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the deviation of B ~ M 1 ' 2 is at most AB ~ (M + AM) 1 / 2 - M 1 / 2 ~ AT" 1 / 2 AM. On 
the other hand, we saw in the last section that p(u) ~ ]xl(M / Bq)\J B 2 — w 2 . Then the 
second integral is suppressed by ln(M / B )B 1 ^ 2 (AB) 3 ^ 2 ~ M -1 / 2 InM, which vanishes 
when M is sent to infinity. After all, in the limit M — > oo we obtain 



H»H*) + 2f-^. (4.7) 

J-oo i U - V > + 1 

This integral equation is solved in the Fourier space [19,36]. Using the techniques 
in the Appendix D of [19], one can derive the following formulas^ 

poo poo 

g 2 /x ± (u)x ± (u 1 )) = 2g 2 dte ±iut e~ t/2 dtfe ±iu ' t 'e- t ' /2 H m (2gt, 2gt'), (4.8) 

Jo Jo 

poo poo 

In (1 - g 2 /x ± {u) x T (u)) = -2g 2 / dte ±iut e~ t/2 / dt'e*™'*' ' e~ v ' /2 K m (2gt, 2gt'), (4.9) 

Jo Jo 

where the integral kernels are expressed in terms of Bessel functions J n (t) by 

ii m (t,t) = , K m [t,t) = — — . (4.1U) 

With these formulas, one immediately obtains the solution in the Fourier space 

i r°° 

a(u) = — / dte ttu a(t), (4.11) 

where 

2 f N l roo 

a(±t) = -4i-j V\ / dt ' e ~ t,/2 (e ±it,w ^H m (2gt, 2gt') + K m (2gt, 2gt') 

S1Iin 2 j=i Jo V 

(4.12) 

for t > 0. 

4.2. The dressing phase 

We are now in a position to compute the dressing phase. The vacuum configuration in 
the bare description corresponds to the empty state in the physical description. The 
physical fundamental excitations are described by adding extra roots w^k in the bare 
description. The bare configuration we studied in the last subsection corresponds to 
the system of N± excitations. 



6 We assume u, u' G R. The branch of logarithm should be chosen appropriately. 
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In the physical description, the scattering phase of the two fundamental excitations 
is simply given by 

U>4.2 + i 



'12 



1 , ^4,1 

- In — — 

i W41 



+ 26»p h y S (w 4i i,M;4 i 2). 



(4.13) 



U>4,2 - % 

In the bare description, the same scattering phase is expressed by the difference of two 
phases 

0i2 = S 12 (w 4 ,i) - Si(w 4 ,i). (4.14) 

Here 612 is the total phase which the first excitation gains when moving around the 
chain in the presence of the second excitation. 5\ is measured in the same way but in 
the absence of the second excitation [32]. 

The total phase is the phase of the transfer matrix eigenvalue and thus can be read 
from the r.h.s. of the central Bethe equations (I2.16p . By substituting (13. 6p (and the 
corresponding relations for ui t k, U2,k), they read 

J N 4 . 2M 

W^k — W4J + I T-r W^k — U4J + I 



Vth 



y*,k 



n 



n 



J7 L k ,J 3= l ,J 

tt 1 ~ 9 2 /ylk x jj yr l _- dVvlk x 2,j 



9 2 /vtk x tj f = \ 1 - 9 2 lvXk x 2j 



M 



1 



X 



9 2 IVi,k x t,j 



M 



3=1 - 9 2 /VA,k x e,j j: 1 
The total phase is then expressed as 

JV 4 . , 2.W 

i 

Sl...N 4 (w*,k 



9 2 /V4,k X 6,j 

9 2 /yt k x li 



(4.15) 



, ^4 



W^k — W±,j + % 



WAk 



+ 

W4 j — I I ' ' W4 k — U4 j — I 

-g 2 /y7 1 x + (u) 1- g 2 /y7 l x~(u) 



- 2M 

7» 



U>4,fc — U 4 j + 1 



+ - 



hi 



1 -g 2 /yt 1 x+(u) 1 -g 2 /yt 1 x (u) 



p ex (u)du. (4.16) 



Therefore (I4.14p gives rise to 



'12 



1 . U>4,i - W40 + i 

- In — ; ; ; 

i W4 1 — W4 2 — i 



2 

+ - 



In 



g 2 /y4 A x + (u) 1 -g 2 /y 41 x (u) 



where (Ti 2 (w), (Ji(u) are given by (I4.lip - fl4.12p with N4 = 2,1, respectively. The first 
term is the bare scattering phase which the first excitation directly feels against the 
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second one. The second term comes from the scattering of the first excitation against 
the stacks whose density deviation 012 — <J\ encodes the back reaction from the sec- 
ond excitation. From the comparison with (14.131) . we see that the second term plays 
precisely the role of the dressing phase. By using the formulas (I4.8l) - (l4.9p again, the 
second term can be expressed concisely in the Fourier space 

20phy S (w 4 ,l,W4,2) 

/oo poo . 

dtj^e* dt'e u ' W4 ' 2 e~ V (jt{2gt,2gt') - £{2gt',2gt)) , (4.18) 

where 

JC(2gt, 2gt') = Ag 2 / dt"K m (2gt, 2gf)- p -—H m (2gt", 2gt') . (4.19) 
Jo e — l 

This precisely agrees with the BHL/BES dressing phase (12.61) . 
5. Discussion 

We have computed the two-body scattering phase of the fundamental excitations over 
the physical vacuum, which precisely agrees with the BHL/BES dressing phase. By 
taking account of the centrally extended psu(2|2) © psu(2|2) symmetry, this suffices to 
determine the whole 256 x 256 components of the S-matrix. From this S-matrix, one 
can construct the Yang equations and derive the complete set of physical Bethe ansatz 
equations, as explained in Section 2. Altogether, our formulation proposes a derivation 
of the asymptotic all-order Bethe ansatz equations with the BHL/BES dressing phase, 
purely based on the symmetry and the integrability. 

In the last section we have considered particular excited states that consist of only 
N4 excitations. For these states the correspondence between the bare description and 
the physical one may be trivial in the sense that each extra bare root represents a 
physical root. However, the correspondence is not so simple in general: When the 
occupation numbers of the physical roots are still equal to that of the extra bare roots, 
the values of the roots could differ. In general, addition of extra Bethe roots at nested 
levels partly breaks the structure of the stacks. A single physical root sometimes 
corresponds to a complex of bare roots and holes. It would be interesting to clarify the 
correspondence. 

For the moment we do not know whether a Yang-Mills operator of finite length can 
be directly realized in the bare description. A possibility is that a physical operator 
of length L ex could be expressed by a state in the chain of length L + L ex while the 
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vacuum is denned in the chain of length L. Of course both L and L + L cx have to be 
sent to infinity, but still the difference would make sense. 

We have determined our vacuum configuration as the simplest consistent solution 
that generates the BHL/BES dressing phase. However, ultimately we wish to derive it 
as the ground state of a certain Hamiltonian. The Hamiltonian has to be expressed in a 
form compatible with the bare description, preferably in terms of the su(2|2) R-matrix. 
It may be derived from the gauge-fixed light-cone Hamiltonian for the Green-Schwarz 
superstring theory in AdS$ x S* 5 [6], which exhibits the invariance under the centrally 
extended psu(2|2) ©psu(2|2) symmetry when the worldsheet is decompactified [7]. 

Our microscopic formulation will be of fundamental use in various directions under 
the latest investigation, for example, the boundary S-matrix [52], the wrapping inter- 
actions [53,54] and the Baxter equations [55,56]. We hope to report the progress in 
these topics elsewhere. 
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the 12th Claude Itzykson Meeting held at the Service de Physique Theorique, CEA 
Saclay and at the Ecole Normale Superieure for their kind hospitality. 

Appendix A. Effective momentum phase 

The centers of stacks u 6 ,k obey the following effective Bethe equations 

m . 
p <*(«6, fc ) = TT M 6,fc ~ u 6j + i 



Acknowledgments 



where 



«6,fc)> 



(A.2) 



1 



2M 




(A.3) 



Note that 




(A.4) 
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where the index j is defined as 



j = j-M-l 



= -M + §, -M +§,..., M-f,M-§. 
Let us evaluate the function in the large M limit: 

2M 



(A.5) 



$ + (u) 



1 Vln l ~ g2 / x ^ x+ ^ 
i ~! 1 — 9 2 l%\ j x + {u) 



1 



M-l/2 



7 £ ln 



y/(j-h) 2 + *9 2 + (J " s) + ( u ) 



A/+1/2 J(J+ |) 2 + V - (J + |) - 2z£+( M ) 



M-l/2 



- V In 

J=-M+l/2 



t? + a 2 - tj + 6 



tf + a 2 - tj - b 



where 



U 



J + 



2l 



a = — , b 



2ix + (u) 



3 M ' M' M 

In the large M limit, one can approximate the sum in (IA.6I) by the integral 



.1/ 



'^i Vt T T^-t + b 

at In , 

i Vt 2 + a 2 - t - b 



(A.6) 



(A.7) 



(A.8) 



This integral can be performed by the change of variable s = \Jt 2 + a 2 — t. In fact, 
j dt\n(Vt 2 + a 2 -t + b) 

Ln 2 \ 

ln(s + 6) 



v / I+^'-l 



Vl+a 2 +l 
2 



s 2 + a 2 



2s 2 



(— - ln(s + 6) - bki(s + b) + y ln(l + - j + s + 6 



=Vl+a 2 -l 



s=Vl+a 2 +l 



(A.9) 



Using this, one obtains 



2i In 



(A + iyl + a 2 ) 2 + 1 
(A - iv / TTa 2 ) 2 + 1 

A 2 + (VTTa 2 + l) 2 
A 2 + (v / TT^ 2 -l) 2 



+Aln 



a 2 , a 4 /A 2 + (v^T^+l) 2 
— In 



A 



a 



VA 2 + (v / T+a^-l) 2 



,(A.10) 
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where 

•i 2x + (u) 2g IK i . 

\ = -ib= — -jr-L a=v?. (A.ll) 

takes the same form with A = 2x~(u)/M. 
Let us consider the case where the coupling constant g is finite. As we take M and 
u sufficiently large, we see that 

2u 

A^— , a^O. (A.12) 
As a result, the phase function reduces to a reasonably simple form 



$(u) = 2M 



r2 



u u . ( M 

2 arctan — + — In H - 

M M V u 



2nM. (A.13) 



In the above computation we implicitly chose the branch of logarithm so that = 
at u — +oo. In the main text we drop the constant —27rM, which corresponds to the 
choice of the branch where $(tt) = at u = 0. 
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